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Abstract 

We study heteroclinic standing waves (dark solitons) in discrete nonlinear Schrodinger equa- 
tions with defocussing nonlinearity. Our main result is a quite elementary existence proof for 
waves with monotone and odd profile, and relies on minimizing an appropriately defined energy 
functional. We also study the continuum limit and the numerical approximation of standing 
waves. 
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1 Introduction 

Discrete nonlinear Schrodinger equations (DNLS) are prototypical models for the dynamics of 
coupled oscillators and have a broad range of application in nonlinear optics and particle physics, 
see [KRB01 , IEJ03|, lKev09al IPor09j . A special focus in the mathematical analysis of such systems 
lies on the investigation of coherent structures (travelling or standing waves, breathers, kinks) as 
these describe the fundamental modes for energy localization and wave propagation in nonlinear 
discrete media. 

This paper concerns heteroclinic standing waves for the one-dimensional DNLS equation 

±Aj - p(A j+1 + Aj_ x - 2Aj) + y'(\Aj\ 2 )Aj = 0. (1) 

Here Aj = Aj(t) denotes the complex amplitude of the oscillator j £ Z at time t and f3 is the 
coupling constant. The potential function ^ describes the energy of a single oscillator system and 
is often assumed to be monomial. In what follows we consider an infinite number of oscillators and 
set either Z = Z (on-site setting) or Z = | + Z (inter-site setting). The difference between both 
settings will be explained below. 

The lattice system ([1]) is an infinite dimensional system of Hamiltonian ODEs (well-posedness 
results are given in [PKM05j IGPlOj ) and possesses the conserved quantities 

H(A) = J>( L4|f ) + \A j+1 - A/ , M(A) = £ \A 3 \ 2 . (2) 
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More precisely, T~L is the Hamiltonian of ([I]), whereas Af gives the power of a state A and stems 
from the gauge symmetry, that is the invariance of ([T]) under the transformation Aj(t) t— > e lip Aj(t) 
with ip e K. 

Standing waves are exact solutions to (pQ) that satisfy = e 1CTt Uj(t) with frequency ffgK and 

real-valued profile u = [uj).^ z G £°°(Z). Standing waves can be regarded as relative equilibria 
with respect to the gauge symmetry and satisfy 

auj = —/3(uj + i + Uj_i - 2uj) + ^'(Uj)uj. (3) 

Heteroclinic waves connect different asymptotic states via 

j— >±oo . . 

uj > u ±oc (4) 

and correspond to so called dark solitons. Other types of standing waves are periodic waves with 
Uj = Uj + N for some N, and bright solitons, which are homoclinic with limj_ 5 .-|- 00 Uj = 0. 
In this paper we aim in establishing the existence of standing waves in the set 

M = {u £ £°°(Z) : u_j = —Uj, — Uqo < Uj < Uj + i < Uqo V j G Z}, 

which consists of all profiles that are odd, increasing, and take values in [— u^, itoo]. In particular, 
we suppose < Uqo = — ii_oo. Notice that we have no = for all on-site waves, whereas inter-site 
waves generically satisfy u_i/ 2 < < ui/ 2 . 

For convex \& one distinguishes between the focussing and the defocussing case, which (in our 
notation) correspond to /3 < and f3 > 0, respectively. Although both cases are linked via the 
staggering transformation Uj (—l) J Uj they describe different physical situations. In the context 
of standing waves it is well established, see |Kev09bj and |Kev09c| . that the most fundamental (that 
means most stable) waves are homoclinic for (3 < but heteroclinic for /3 > 0, respectively. In 
what follows we solely consider the defocussing case (3 > and construct heteroclinic single-pulse 
waves by minimizing an energy functional. The analogues in the focussing case are homoclinic 
single-pulse waves, which can be constructed by constrained maximization, see |Herl0j . 

Several methods have been used to prove the existence of standing waves in DNLS equations. 
Continuation methods were introduced by MacKay and Aubry |MA94|, IAub97j and have been 
proven powerful for both theoretical considerations and numerical computations, see for instance 
|KL09j . Continuation method rely on the observation that © can be solved explicitly in both the 
anti-continuum limit /3 — > and the continuum limit /3 — > oo. During the last years, however, 
there has been a growing interest in other approaches to the existence problem for standing waves. 
We refer to |PR05| . which exploits spatial dynamics and centre manifold reduction, and to the 
variational methods in jPZOH IPR081 IZP09|, IZL09] , which rely on critical point techniques (linking 
theorems, Nehari manifold). 

We now summarize the main idea in our variational existence proof. At first we notice that Q 
and couple the frequency a to the asymptotic states via 

* = *'(ul). 

Secondly, we introduce the function 

F{n) = *(n 2 ) - *(t&) - {V 2 ~ 
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and rewrite the standing wave equation as 

F'( Uj ) - 2/3(u j+1 + uj-i - 2 Uj ) = 0. 
The key observation is that each standing wave is a critical point of the energy functional £ with 
S (u) = T{u) + PV{ U \ F(u) = Y,F{u j ), V{u) = Y J W 3+ i-u j ) 2 . (5) 

3 j 

The energy £ is naturally related to the conserved quantities from (J2J). In fact, on a formal level 
we find 

£{u) = %{u) - <jN{u) - U{u) + o~j\[{u) + pD{u) 

where u is an arbitrary reference profile with = for all j, but notice that the M- and %-terms 
are infinite if u satisfies (j3J). £(u), however, is well defined as long as u approaches the asymptotic 
states sufficiently fast. In fact, using F^u^,) = F'(±Uoo) = we find that both J~{u) and T>{u) 
are finite if the sequence j i-> Uj — UooSgnj belongs to £ 2 {Z). 

Our strategy for proving the existence of standing waves is to show that £ attains its minimum 
on M by using the direct method from the calculus of variations. Afterwards we show that each 
minimizer satisfies ^ as it is strictly increasing. Our existence result for energy minimizing waves 
can be summarized as follows. 

Theorem 1. Let > be given and let ^> be twice continuously differ entiable on [0, x^] with 
= u^q. Moreover, suppose that 

^"(xoo) > and W(x) > $(x) + ^'{x 00 ){x - x^) for all < x < Xoo- (6) 

Then, for each (3 > the functional £ attains its minimum on Ai. Each minimizer u £ M is 
strictly increasing, converges exponentially to ±Uoo as j — > ±oo, and solves the standing wave 
equation ([3|) with frequency a = ^>'(xoo). 

We proceed with some remarks concerning the assumptions and assertions of Theorem [TJ 

1. Theorem [1] holds in both the on-site and inter-site setting. Moreover, via Uj ~^ —Uj it also 
provides also the existence of standing waves with decreasing profile. 

2. Suppose that is strictly convex on [0, oo) with ^'(0) = 0. Then we have ^f(x) — ^(y) > 
*$>'(y)(x — y) for all x, y > and Theorem [1] guarantees the existence of standing waves with 
arbitrary u^. Furthermore, there exist standing waves for non-convex \E r , see the examples 
in g3j 

3. Assumption §§§2 precisely means that F is positive on the interval (— Uoo, Uoo). This condition 
implies F'^Uoo) > and is sharp in the following sense. Suppose there exists < < 
such that F(u*) < 0. Then the energy £ is unbounded from below and global minimizers of 
£ can therefore not exist. See also the discussion in g3l 

4. Condition ([6])i is equivalent to -F"(uoo) > and guarantees that each standing wave ap- 
proaches its asymptotic states exponentially. For F"{u 00 ) = one can still prove the exis- 
tence of energy minimizing waves (using some limit procedure) , but since £ has no minimizer 
for F"(u OQ ) < these waves are non-generic and have an algebraic tail. 

5. The standing waves from Theorem [1] are so called single-pulse dark solitons. From continu- 
ation results we know that there also exists an infinite number of multi-pulse dark solitons 
corresponding to non-monotone solutions to and see for instance |Kev09c] . However, 
multi-pulse waves are expected to be unstable and we conjecture that they correspond to 
genuine saddle points of £. A detailed investigation of this issue is left for future research. 
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The proof of Theorem Q] is given in E j2.ll and Q2.21 In £12.31 we establish an approximation result 
which allows to compute standing waves on finite index sets and provides the base for the numerical 
simulations in ^3 Moreover, in £12.41 we study the continuum limit of standing waves with fixed 
asymptotic states. To this end we introduce a scaling small parameter e > and scale §5§ by 
j ej and /3 ~» /3/e 2 . Then we proof that the standing waves converge as e — > to a heteroclinic 
solution of 

2Pu"(0 = -F'(u(0), .lim u(0=±«oo- (7) 

2 Existence and properties of standing waves 

In order to prove Theorem [1] we first observe that ([3]) is invariant under the scaling 

Uj i y rjUj, a i y to~, f3 h- )• r/3, ^(x) i-» Tri~ 2 ^{j] 2 x) + 8 
with arbitrary 77, r, <5 G R. Prom now on we therefore assume that 

a = *(1) = = 1, u ± oo==bl, F(r?) = ^(r ? 2 ) -r/ 2 , 

and study the normalized standing wave equation 

f rf («,-)-2 J 8(«,- +1 + Ui_i-2u i ) = 0, Iim«j-=±1. (8) 

j— >±oo 

2.1 Standing waves as minimizers of £ 

Our assumptions on ^ ensure that £ has nice properties on Ai. Recall that Ai is compact with 
respect to the weak* topology in £°°(Z). 

Lemma 2. £ is non-negative and weakly* lower semi- continuous on Ai. Moreover, £{u) < 00 
implies lim Uj = ±1 and 

lim + to) - £(u)) = Vj Qlu), (9) 

/or all v G ^(Z), w/iere Q(u) G r°(Z) tuitt 

= "S>'(uj)uj - Uj - (3(uj + i + uj-i - 2uj) 
is the Gateaux derivative of £ in u. 

Proof. Let u G Ai be given. Then we have £ (u) > J-(u) > due to the assumption on i* 1 , and 
the monotonicity of u implies lim J _ > .-i- 00 Uj = ±Uoo for some Uoo G [0, 1]. Assuming Uoo < 1 we find 
F{ u j) > \F( u oo) > for almost all j, and hence £{u) > .F(ii) = 00. Moreover, ([9]) follows from 
a direct computation. Finally, we consider a sequence (u n ) neN C .M, with u n = (u n j)j e z ^ or an 
n, that converges weakly* in £°°(Z) to some limit u G «M. Then we have u n — >■ Uj for all j, and 
Fatou's Lemma gives £(u) < liminf n _ 5 . 00 £(u n ). □ 

We are now able to prove our main result. 

Theorem 3. £ attains its minimum on Ai. Moreover, each minimizer is strictly increasing and 
satisfies the standing wave equation ([8|). 
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Proof. The existence of minimizers follows from the weak* compactness of Ai and the lower semi- 
continuity of £. We now show by contradiction that each minimizer u G Ai is strictly increasing. 
Suppose at first that there exists some jo G Z such that 

uj -i < 1 and Uj = 1 for all j G Z with j > jo, 

and define v £ ^(Z) hy Vj = — 1 and Vj = for all j ^ jo- Then we have 

G(u) jo = tf'(l)l -1-/3(1 + %0 _! - 2) = /3(1 - u j0 -!) > 0. 

and ([9]) yields £(u + tt>) < £{u) for all sufficiently small t > 0, which is the desired contradiction. 
Consequently, we have Uj < 1 for all j £ Z and —1 < Uj follows similarly. Secondly, suppose there 
exist a constant plateau between —1 and +1, that means 

u ji-i < u h = = ••• = u h-i = u h < u j2+i 

for some indices j\ < j2, and define v G ^(Z) by Vj 1 = —1, fj 2 = +1, and vj = for all 
j G Z \ {ji, j 2 }. This implies 

G{u) h = ^'(u 2 h )u h - u h + P{u h - u h -{) > $'(uj 2 )u h - u j2 + /3(u j2 - u j2+1 ) = G{u) h , 

and (J9j) gives £{u + tv) < £(u), which is again a contradiction. Finally, since u is strictly increasing 
we can perturb the itj-'s independently from each other, so ^ yields G(u)j = for all j G Z. □ 

Recall that Theorem [3] provides the existence of standing waves in both the on-site and off-site 
setting. 

2.2 Exponential tails for standing waves 

We next show that standing waves converge exponentially to the asymptotic states. Heuristically, 
the decay rate A is determined by linearizing ([8]) in the asymptotic states, i.e., A is the unique 
positive solution to 

4/3(coshA- 1) = F"{1) = 4tf"(l). (10) 

From this we conclude that \&"(1) > is truly necessary for standing waves to have exponential 
tails. 

Lemma 4. Let u G Ai be any strict monotone solution to (J8j) and choose X, A such that < A < 
A < A with A as in (|1U[) . Then there exist constants c and c such that 

cexp(-A|j|) < |sgnj-nj| < cexp(-A|j|) 

for all j G Z. 

Proof. By virtue of U-j = —Uj it is sufficient to consider positive j. For 1 < j G Z we define Wj 
and Kj by Wj = 1 — Uj and Kj = Wj/wj-\. By u G M. and ([8|) we have 

< Wj < Wj-i < 1, 5jWj = Wj + i + Wj-i — 2wj, 5j = — J , 

2p 

where £j denotes some intermediate value in [uj, 1], and hence 

< Kj < 1, Kj + l = 2 + Sj . 
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Moreover, using Wj — > as j — > oo we find 



- j-^oo . F"(l) j^oo v _ 2 + 5 oc - y^y^ + dpo 

Oj > Oqo — , Kj > Kqo — ^ ) 

and a direct computation reveals that A = — In (k^) > satisfies (fTUj) and A < A < A. The desired 
result now follows immediately. □ 



We mention that (|10p implies the expansions 

F"(l) 

2(3 ' 



A = -ln(5~ 1 ) + 2<T 1 (l + 0(<r 1 )), A = v^(l + 0(<5)), <5 



In particular, we have A — > oo in the anti-continuum limit (3 — > but A — > in the continuum limit 
P -> oo. 

2.3 Ritz approximation of standing waves 

It is reasonable, and useful for numerical simulations, to approximate the standing waves from 
Theorem [3] by minimizing £ on the finite-dimensional set 

M N = {u £ M : Uj = sgnj V j G Z with \j\ > N}, N G N. (11) 

Notice that £ attains its minimum on .Mat as .Mat is closed under weak* convergence in £°°(Z). 

Lemma 5. We have 

mmt\M N > mmc\M- (12) 

Moreover, each sequence (un) N( - n with un G Mn and £(un) = mm£\M N for all N, has a 
subsequence that converges pointwise to some u with £(u) = m.m.£\_M. 

Proof. Let u G M. be a minimizer of £ , and define the sequence (itjv)jv by u iV,j — % for |j| < N 
and UN,j = sgnj for |j| > TV. Then we have 

Hun) = Yl ^ 
b'[<Jv 



and 



V(u N ) = {uj+i-Ujf + 2{l- u N f <V{u) + 2(1- u N f. 

\i\<N-i 

Taking the limit N — > oo gives 

lim min<f < lim £(un) < £( u ) = min5|x- 

N^oo N— >oo 

This implies ([12]) because we have M.n± cA4n 2 CM and hence minf^ < minfl^^ < min£ \m n 
for all N\ < N2. Finally, suppose that (un) n satisfies £(un) = m.m£\M N . By weak* compactness 
we can extract a (not relabelled) subsequence that converges weakly* in £°°(Z) to some limit 
uGjM. Lemma [2] combined with (|12p gives 

min£|x < 5 (it) < lim £(un) = lim min£|;\4^ = min£|_A/f , 
so it is in fact a minimizer of £. □ 



6 



2.4 Continuum limit of standing waves 

In this section we characterize the continuum limit of discrete standing waves. To this end we 
consider profile functions u = it(£) G M, where £ G R is continuous variable and Ai is given by 

M = {u G L°°(R) : -1 < n(£i) < = -«(-&) < 1 V & < & € R} . 

The set M is a naturally related to the formal limit problem fl7|) as standard arguments for planar 
Hamiltonian ODEs imply the following result. 

Remark 6. There exists a unique u G A4 with F'(u) = 2(3u" . 

We now introduce a small lattice parameter e > and consider the functional £ £ in M defined 

by 

£ £ (u)=F(u)+(3V £ (u), T(u) = jF(u(0)dC, V £ {u) = J M^lAz^ J ^, 

R R 

The connection with the discrete setting becomes obvious when introducing 

M £ = {u G M : u(ej) = u{ej+et) V j G Z, £ G ±) G R} . 

Each it G A4 £ is piecewise constant and can be identified with a increasing and odd profile from 
£°°(Z). In particular, we have 

F{u) = e v e(u) = ^ + e ) " U ^)Y 

jcz jez 

The result from §2.11 imply that £ £ attains its minimum on M £ . 

Lemma 7. For each e > there exist a minimizer u £ G M s for £ £ \m e such that 

J F'MOMOV = V / M yK + e) + y( 4 + e) " iM0 di (13) 

R R 

ZioWs /or all test functions ip G Cg°(R). Moreover, there exists a constant C such that £ £ {u £ ) < C 
for allO < e < 1. 

Proof. The existence of u £ follows from Theorem (with /3/e 2 instead of and (|13p is the weak 
formulation of the discrete standing wave equation, see ([8j). We now consider v £ G Ai £ with 



— 1 for all j G Z with ej < —1, 
+1 for all j G Z with ej > +1, 
ej for all j G Z with e | j| < 1. 



A direct computation shows 

i 

F(v £ ) = Y, ^ F (ej) = [ F( v ) dr, + O(e), V £ (v £ ) = ^ e = 2 + 0(e). 

e\j\<i i e\j\<l 

The uniform bounds for £ £ (u £ ) now follows due to £ £ {u £ ) < £ £ (v £ ). □ 
We are now able to pass to the continuum limit e —¥ 0. 



Corollary 8. We have u £ — > u as e — > in Lj^ c (R), where u £ and u are defined in Lemma^7\ an d 
Remark® respectively. 

Proof. Let (e n ) n C M be any sequence with e n — > 0, and let u n = u £n G M En . By weak* 
compactness in L°°(1R) we can extract a (not relabelled) subsequence that converges weakly* to 
some limit Uqo G Ai. Lemma [7] guarantees that the sequence (v n ) n with 

u n (£ + e n ) - u n {(,) 

is bounded in L 2 (R) , and passing to a further subsequence we can assume that v n converges weakly 
in L 2 (M) to some limit Voq. We then have 

= J v n (0v(0 d£ + J u n {i f {i) ~^~ £n) d£ V V9 G C °°(M), 

and passing to the limit n — > oo we conclude that Voo is the weak derivative of Uoq. In particular, 
Uqo is continuous and satisfies 



5 ? 
UooiO = / v OD (()d(= lim / Vn (()d(= lim u n (£), 

J n— ¥oo J n— ¥oo 




where we used that ii n (£) = ^ v n {C) dC, + 0(e n ) holds by construction. Therefore, u n converges to 
Uoo in LJ^ C (IR), so f)13|) implies that Uoq and u solves the same ODE. Moreover, Lemma [7] combined 
with Fatou's Lemma provides 

Flu) < liminf F{u n ) < oo, 

n— >oo 

which implies Urg* — > ±1 as ^ — > ±oo and hence = u. We have now shown that each sequence 
(u £n ) n has at least a subsequence that converges to u. The uniqueness of u now implies the desired 
result. □ 



3 Numerical simulations 

We illustrate the analytical results from the previous section by numerical simulations of discrete 
waves. To this end we implemented the explicit Euler-scheme for the gradient flow of £ in the Ritz 
set Mn from (jlip . The corresponding iteration scheme with flow time r > reads 

u G Mn !-> 1{u) G M n, 

with 

X{u) j = Uj - r{F'(u j ) - 2p(u j+1 + Uj-x - 2 Uj )) V \j\ < N. 

The iteration scheme works very well in numerical simulations and converges to a fixed point of 
I provided that r is sufficiently small. It is clear, at least for large N, that each fixed point of 
X satisfies the standing wave equation up to high order, but it may happen that this fixed point 
approximates a local minimizer of £. 

To ensure that the scheme converges to a global minimizer we need a good guess for the initial 
profile. In our simulations we always start with the shock profile Uj = sgnj for all j, which is the 
global minimizer for f3 = 0. 
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Figure 1: On-site and inter-site waves for /3 = 0.25 after 150 steps with r = 0.1 and N = 6. The 
upper left picture shows the graphs of tp (Black) and ^ (Gray), see (Q2 
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Figure 2: On-site waves for j3 = 1 and [3 = 5 after 600 steps with r = 0.01 and N = 12. 



Typical numerical solutions for convex and non-convex \T/ are shown in Figure [T] and [21 respec- 
tively. In the upper left picture we also plotted the graph of 

V>(r?) = *'{v 2 H (14) 

which is the nonlinearity in the standing wave equation Notice that the first example is 

prototypical for power laws V ; ( 7 ?) = T r+^ r / 1+d with d > 1. 

Figures [3] illustrates what happens if ^ violates Q. In this example, F is not positive in (—1, 1) 
but there exists < r/* < 1 such that > ^(t?*) = F(— 77* ) = min-i^^i F(rj). Consequently, £\m 
is unbounded from below and global minimizers u £ Ai cannot exit. In fact, setting Uj = i]* sgnj 
for all \j\ < N and Uj = sgnj for \j\ > N we find £ (u) = NF(7]*) + O(l) — > —00 as N — > 00. We 
now interpret the numerical results in Figure El For small (5 (and large N) the iteration scheme 
converges to a fixed point of I that satisfies the standing wave equation up to high order and has 
a sharp transition from —1 to +1. We therefore conjecture that standing waves with it±oo = ±1 
still exist and correspond to local minimizers of £\m- F° r larger values of f3, however, the solutions 
un £ M-N exhibit plateaus at heights ±77* and converge for N — > 00 to some limit profile u with 
asymptotic states ii±oo = Notice that u solves the standing wave equation ([8]) as F'(rj*) = 
implies ^'(t?*) = 1. 
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Figure 3: Off-site waves for /3 = .5 and f3 = 2 after 1000 steps with = 40 and r = 0.05. 
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Figure 4: Continuum limit of inter-site waves with (3 = 1: Discrete waves (Black) from A4 e for 
several values of e along with the solution (Gray) of the limit problem ([7J). 



Finally, Figure[5]concerns the continuum limit of standing waves, see £ )2,4I It shows the piecewise 
constants minimizers u £ € M e for different values of e in the interval |£| < 6, and illustrates that 
u £ converges for e — > to the unique solution of ([7D . 
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